Terzan 5: an alternative interpretation for the split horizontal branch 
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ABSTRACT 

We consider the horizontal branch (HB) of the Globular Cluster Terzan 5, recently shown 
to be split into two parts, the fainter one (5Mk ~ 0.3mag) having a lower metallicity than 
the more luminous. Both features show that it contains at least two stellar populations. The 
separation in magnitude has been ascribed to an age difference of ~6 Gyr and interpreted 
as the result of an atypical evolutionary history for this cluster. We show that the observed 
HB morphology is also consistent with a model in which the bright HB is composed of sec- 
ond generation stars that are metal enriched and with a helium mass fraction larger (by 
5Y ~0.07) than that of first generation stars populating the fainter part of the HB. Terzan 
5 would therefore be anomalous, compared to most "normal" clusters hosting multiple pop- 
ulations, only because its second generation is strongly contaminated by supernova ejecta; 
the previously proposed prolonged period of star formation, however, is not required. The 
iron enrichment of the bright HB can be ascribed either to contamination from Type la 
supernova ejecta of the low-iron, helium rich, ejecta of the massive asympotic giant branch 
stars of the cluster, or to its mixing with gas, accreting on the cluster from the environment, 
that has been subject to fast metal enrichment due to its proximity with the galactic bulge. 
The model here proposed requires only a small age difference, of ~100 Myr. 

Subject headings: Galaxy: bulge — (Galaxy:) globular clusters: individual (Terzan 5) - 
stars: horizontal-branch — stars: AGB and post-AGB — supernovae: general 



1. Introduction 



Recent years have witnessed exciting developments both in the observations and theoretical mod- 
elling of the abundance star to star variations within most of the well studied Globular Clusters (GCs). 
In most GCs star to star abundance variations are limited to the light elements that are susceptible to 
abundance changes from proton- capture reactions , such as the pp, CN . ON, NeNa, and MgAl cycles 
(for t he most recent spec troscopic survey, see, e.g. Carretta et al.J 2009al Jbl). but a few clusters, such as 



M22 ( jMarino et al.l 120091 ). or perhaps NGC 1851 (Han et al.ll2009l ) are now known to exhibit variations 
in heavier elements (see also ICarretta et al.ll2009a). arid, m ore than the others, in wCen the heavy ele- 
ments spreads (e.g., among others iNorris fc Da Costalll995l ). and the HR diagram morphologies clearly 
show that we are dealing with several st ellar generations, enriched by the supernova (SN) ejecta (e.g. 



Sollima et al.ll2005l ; IVillanova et al.ll2007l ). In addition, the cluster M 54, immersed in the nucleus of the 
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Sagittarius dwarf galaxy, presently disrupting in our Galaxy (llbata et al.lll994t iBellazzini et al.l 120081 ) 
has been recently found to show a metallicity spread similar to u Cen (jCarretta et al.ll2010l ) 



Concer ning the spread in l ight elements, their observation at the turnoff and among the subgiant 



stars (e.g., iGratton et al.l 120011 ) showed that these anomalies must be attributed to self-enrichment 



occurring at the first stages of the life of the cluster. The most peculiar finding of the latest years is the 
presence of a very helium rich population in the most massive clusters: this is revealed by the p resence 
of multiple main sequences in ui Cen and NG C 2808, indicating a helium content Y=0.38-0.40 (INorris 
2004 ; D'Antona et al. 2005 ; Piotto et al.l 2007 ). and by the extreme morphology of the horizontal branch 
(HB) of the two massive clusters NGC 6388 and NGC 6441. In these latter cluste rs, a red clump is 
expected as HB, due to their large metallicity ([Fe/H]~-0.4, see Carretta et al.l 2009c ). On the contrary, 
their HB is extended tow ards the blue, and t he RR Lyr variables have so long per iods that they must 
be highly overluminous. Caloi fc D'Antona ( 20071 ) and D'Antona &; Caloi fl2008t ) show that the HB 
morphology and RR Lyr's periods of these two clusters may be explained if a large fraction of the HB 
stars have Y £ 0.35. The presence of much more moderate helium spreads is probable in most of the 
other smaller clusters (e.g. iD'Antona fc Caloill2008l ). The quasi-constancy of heavy metals in most GCs 
leads to hypothesize that the abundance variations must be due to very peculiar chemical evolution, not 
or scarcely affected by SN ejecta, but involving formation of a "second generation" (SG) of stars from 
matter processed into the "f irst generation" (FG) stars. On the other hand , the numerical consistency 
of the SG is so high (>50% ID'Antona &; Caloil 120081 ; ICarretta et al.ll2009al ) that any formation model 
must include the hypothesis that the mass contained in FG stars — that contribute to this second phase 
of star formation — is initially much larger than the mass present today in the cluster. Models for the 
formation of these multiple generations are still in their infancy. We can divide them into two main 
categories: the models in which clusters are born inside dwarf galaxies, so that t he polluting matter 
on the forming GC comes from a much larger environment (e.g. iBekki et al.l 120071 ) . and the models in 
which there is an initial cluster 10-20 times more massive than todays'. In the latter case, the SG forms 
from the ejecta of the FG stars mainly in the centra l cluster parts, and t he first dynamical phases of 



evolution lead to a preferential loss of the FG stars ( D'Ercole et al.l 2008 ). In these models, it is very 



difficult to accomodate large age differences between the first and second generation stars. 

Consequently, the recent observations of the color magnitude diagram features and chemistry of 



Ferraro et al. 



Terzan 5 may constitute a benchmark in our understanding of GC formation. In fact 
(120091 ) show that the cluster HB stars are divided into two clumps separated by 5Mk ~ 0.3mag, and 
that the more luminous stars have a much larger iron content ([Fe/H]~+0.3±0.1) with respect to the 
lower HB ([Fe/H]~-0.2±0.1). This result shows that the evolution of this cluster is atypical, and that 
matter forming the SG stars (populating the u pper HB clump), has been affected by SN contamination, 
as it occurred in u Cen. On the other hand, iFerraro et al. I (l2009h . comparing the HB data to stellar 
isochrones of the correct metallicity, conclude that the SG must be ~6 Gyr younger than the FG. This 
huge age difference is very difficult to be understood in any formation fram ework, and this would be 
the first evidence for such a young age among bulge stars and clusters (e.g. iFeltzing fc Gilmord 12000 
Origlia et al.1l2008h . 



We re-examine the problem and show that the HB morphology can also be explained by two coeval 
populations having an helium difference of <5Y~ 0.07, thus reaching a value not as extreme as in the 
cases quoted above, so its formation does not present particular problems (Sect. SJ). At the super-solar 
metallicity of the SG of Terzan 5 the possible helium enrichment in the SG does not produce a blue 
extension of the HB. In Sect. H]we discuss some possibilities for the chemical evolution of the cluster. We 
finally remark that the different space distribution of the two populations might imply a mass difference 
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between them, and that further observations and dynamical modelling may allow to choose between 
models based on age- or helium-differences. 



2. The stellar models 



We computed evolutionary tracks, isochrones and HB tracks for a metallicity Z=0.01 (having 
[Fe/H]~-0.2 for a solar-scaled mixture) and helium content Y=0.26, representing the FG of Terzan 5 
and its lower HB clump, and models and isochrones for Z= 0.03 (|Fe/H]~+0.3), Y=0.29 and Y=0.40. 
The st andard inputs of our ev olutionary code ATON are used (jVentura et al.ll2009l ). We adopt the opaci- 
ties bv lFerguson et al.l (120051) at tempera tures lower than 10000 K and the OPAL opacities in the version 
documented by Iglesias fc Rogers! (19961). The mixture adopted is solar-scaled and follows the element 
distribution by lGrevesse fc Sauval (jl998 ). For comparison, we also computed models for iron content 
[Fe/H l~-0.2 and [a/Fe]=0.4. Electron conduction opacities follow the treatment by iPotekhin et al. 
(119991 ) . and are harmonically added to the radiative opacities. The HB models are constructed by as 



suming the core helium mass derived from the red giant evolution of masses evolving in the range of ages 
10-13Gyr for the given Z and Y. The values assumed are then M C =0.4793M Q for Z=0.01 and Y=0.26; 
M C =O.464M for Z=0.03, Y=0.29 and M r =0.448M m for Z=0.03 and Y=0.40. From the isochrones we 



derive the mass at the tip of the giant branch for each age and chemical composi tion. These values 



are us ed to build up synthetic models of HB, following the procedure described by iD'Antona fc Caloi 
(120081 ) . The mass lost along the red giant branch and its dispersion are fixed for each given age, in order 
to derive the distribution of HB masses. A further random extraction of the age within the HB lifetime 
allows to populate the synthetic HB. We fix the FG population at the chemistry Y=0.26, Z=0.01, and 
the SG at Z=0.03. The helium content of the SG was allowed to vary between 0.29 and 0.40, but 
we assume the same mass loss for both populations. T he theoretical value s of luminosity and T c g are 
converted into the Johnson-Bessell system by means of iBessell et al.l (119981 ). 



3. The HB morpology of Terzan 5: two possible interpretations 

We start by discussing in this section an example that just serves to illustrate how an increase in 
the helium abundance can explain the difference in K magnitude between the two clump s in Terzan 5. 



Figure 1 shows a representation of the results that allows to graphically visualize both iFerraro et al. 



(120091 ) conclusion and a different interpretation of the data proposed here. In the right panel of Fig.l we 
shows the HB mass versus age relations and on the left panel we plot the HB mass versus the absolute 
magnitude of the ZAHB. Let us try to understand why we need ~ 6Gyr of age difference, if we 
assume (Z=0.01, Y= 0.26) as composition of the faint clump and (Z=0.03, Y=0.29) for the brighter one 
(iFerraro et al.ll2009h . Assume that the age of the faint clump is 12.5Gyr. The evolving mass on the 



HB is represented by the point labelled A (age=12.5Gyr and M=O.646M ). This mass on the ZAHB 
has M.K-— 1-16 (point B). Now we know that the other clump is ~0.3mag brighter, so that we shift 
to M^;=-1.46 (point C). This magnitude needs a HB mass of O.9M (point D), if we are dealing with 



1 We plot the mass at the tip of the red giant branch, that is a decreasing function of age, minus 0.3Mq, in order to 
provide an example of what may be the masses populating the HB. 
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Fig. 1. — Right panel: the curves with triangles represent the mass versus age relations for the three 
different compositions (Y=0.29, Z=0.03: top) (Y=0.26, Z=0.01: middle) (Y=0.40, Z=0.03 bottom). 
The mass is the red giant mass evolving at each age, minus 0.3M Q , that represents an average mass 
lost on the red giant branch. Thus the relations represent the average HB mass as function of age. Left 
panel: the curves with squares represent the K-magnitude as a function of the mass in the ZAHBs of 
(Y=0.26, Z=0.01: top); (Y=0.29, Z=0.03, middle); (Y=0.40, Z=0.03: bottom). The letters label how 
we can have two populations differing by 5Mk ~ 0.3mag by increasing Y and Z (from A to G, dashed 
path) or by increasing Z and decreasing the age (from A to E, dash-dotted path) as explained in the 
text. 



models of Y=0.29, Z=0.03. The age corresponding to this evolving mass H is t hen 7Gyr (po i nt E) . 
Thus the brighter clump is at least 5.5Gyr younger than the fainter clump, as Ferraro et al. ( 20091 ) 



find. Let us suppose now that the two populations are practically coeval, and that the brighter clump 
corresponds to a population having Z=0.03 and Y=0.40. In this case, we must have an evolving mass 
of O.528M (point F) that has M#=-1.58 on the ZAHB (point G). In this case, the gap in Mk between 
the two populations is the magnitude difference between points B and G ( — 1.16 + 1.58 =0.42mag). 
This example shows that a helium increase can easily explain the magnitude difference between the two 
clumps, without requiring an age difference. 

Among the many simulations of the HB, produced to understand the role of the different parameters, 
Figure [2] shows the simplest one that reproduces the gap of Mjf ~ 0.3mag between the two populations. 
We assume Y=0.33 for the bright clump. All stars have age llGyr, mass loss along the red giant branch 
<5M=O.28M , with dispersion a=O.O25M . Both the color and luminosity difference of the clumps are 
reasonably reproduced. The choice of models, age and helium enhancement is not unique. The luminous 
clump can also be composed of stars with Y varying in the range Y=0.32-0.34, born from ejecta with 



2 This is the age, if the mass loss does not increase with metallicity, otherwhise it must be considered an upper limit 
to the age 
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Fig. 2. — The dots represent the simulation of the HB in Terzan 5 by assuming a population having 
Y=0.26 and Z=0.01 (fainter clump) and a brighter HB with Y=0.33 and Z=0.03 (upper clump). The 
three lines with squares represent the ZAHBs of our models: Z=0.01 at bottom left, Z=0.03, Y=0.29 
at bottom right, and Z=0.03, Y=0.40 (upper line). 



different helium. If 5Y<0.07 between the two populations, however, the gap can not be reproduced. 
We have assumed for both the FG and SG a solar scaled composition, a reasonable choice for the SG, if 
its higher metallicity is due to Type la SN contamination, that pollute mostly with iron and bring the 
composition towards the solar-scaled abundances. If the FG is instead a-enhanced, and we assume for 
it [Fe/H] =— 0.2 and [a/Fe]=0.4, our models show that the ZAHB shifts to redder J-K, by ~0.03mag, 
and the color difference of the two clumps is reduced. A small reduction of the iron content, within the 
range allowed by the measurements errors, would reproduce again the color fit. 

The iron content of the upper clump may be not unique for all stars, if it is a result of non- 
homogeneous contamination of the matter forming the SG stars. In this case, a proper interpolation 
between HBs of different [Fe/H] must be taken into account, but a similar result will be obtained. 

Notice that both the FG and the SG clumps "stay in the red" , while the other high metallicity 



clusters, NGC 6388 ([Fe/H] — 0.40) and NGC 6441 ([Fe/H] — 0.33) flCarretta et al.ll2009d ). have HBs 
extended towards large T e g- thanks to the much larger Y of their SG. In fact, the iron content of Terzan 5 
is much larger, and the H-burning shell maintains a giant structure and a large radius for HB stars 
even of relatively small mass. Our models show that we need 0.53M Q to depart towards larger T e g for 
Y=0.29 and 0.51M Q for Y=0.40. Notice that a HB mass of O.56M occupies the RR Lyr pulsation 
strip for the FG c hemistry. A single RR Lyr and a few blue HB stars have been found in the cluster 



( jCohn et al.ll2002l ): these may well represent the tail of the mass distribution of HB stars populating 



either of the red clumps, and it would be interesting to know which one. 
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4. The chemical evolution of Terzan 5: type II or type la SN enrichment? 



We have shown that the HB of Terzan 5 may be explained either by two populations with an age 
difference of 6 Gyr or with two populations that are approximately coeval (<5(age)~100 Myr) but have 
a different helium content (5Y~0.07). The latter interpretation would put the formation of Terzan 5 
with in the theoretical fra mework suggested for the formation of most GCs with multiple populations 



e.g. D'Ercole et al. 2008[ ) . However, while in most clusters SG stars have iron abundance similar to FG 



stars, the iron enhancement in the second generation population of Terzan 5 introduces a new ingredient 
in that scenario and requires the identification of the source of metal rich gas. 



We tentatively iden tify the helium rich gas with the massive-AGB and super-AGB ejecta ( [Ventura. D'Anton 
20021 ; iPumo et al.ll2008f). and comp ute the evolution of these stars for the FG composition (Z=0.01), fol- 



lowing l^^ur^&^^Antona ( 2009 ). The smallest star igniting carbon in conditions of semi-degeneracy 
(super-AGB evolution) is the 7.5M , evolving at 50Myr. The AGB ejecta have helium abundance be- 
tween Y=0.36 (7M ) and Y=O.32M (5M ), so masses down to ~5M (evolving at ~100Myr) can fit 
the SG requirements. We have ~100Myr to pollute the AGB ejecta with iron rich material. A change in 
the iron conte nt from |Fe/H]~-0. 2 ±0.1 to [Fe/H]~+0.3±0.1, assuming Z =O.O18, an iron mass fraction 
f=0.074 from iGrevesse fc Sauvall ( 119981 ). and solar-scaled compositions, means an increase in the iron 
mass fraction by 1.8 xlO -3 in all the SG mass. Assuming that this mass is half of the total mass of the 
cluster today, say 2.5xlO 5 M , the requirement is 456M of iron. We can think of two ways of forming 
an SG both helium and metal enriched. 

1) The helium enriched ejecta are directly polluted by iron produced by superno vae belonging to 



the cluster itself. If the source of the helium enrichment are the massive stars, as in Decressin et al. 



( 120071 ) model, the sou rce of the metal enrichment must be the SNII. Each SNII produces only ~O.O7M 
of iron ( iHamuyl 120031 ). so we would need 6500 SNII explosions to reach the required iron content, an 
occurrence that would very likely destroy the cluster, and will not be considered any long er. If the source 



of hel ium enrichment are the massive AGBs, like in the dynamical model proposed by iD'Ercole et al. 
( 120081 ). for Terzan 5 we must further require that the second stage of star formation is not halted by 
the injection of energy from the SN explosions, that generally is regarded as the cause of the end of 
the SG formation epoch. Below 7.5M , as remarked above, stars evolve into carbon oxygen white 
dwarfs, so that, some time after, SNIa can begin exploding in the cluster. Each SNIa injects in the 
intracluster medium ~O.6M of iron, so we need 760 SNIa explosions. A rate of 1 SNIa every 50000yr 
can provide this iron in 38Myr, when the AGB ejecta are still as helium rich as needed. If the SG will 
result to be homogeneous in iron, the model requires that the helium rich AGB gas accumulates, and 
that the SNIa explode polluting the AGB gas, but not destroying its accumulation. The above back of 
the envelope computation shows that the requirements for the increase of the iron content may be in 
principle fulfilled, although with many caveats. In this scenario we expect that in the SG [a/Fe]~0. As 
for the other elemental abundances, this model can be easily falsified, as it predicts that the SG should 
not show the enhanced [Na/Fe] signature typical of GCs. In fact, the AGB gas may well be sodium 
enhanced with respect to its iron content [Fe/H]~ —0.2, but the iron rich mixture forming the SG will 
necessarily have a low [Na/Fe]. 

2) The SG matter consists of the low-iron gas, ejected from the cluster AGBs, mixed with the 
environment gas accreting on the cluster from the neighborhood regions. Due the proximity to the 
galactic bulge, this gas has been subject to fast metal enrichment, and its iron may be largely superso- 
lar even at a very early epoch, depending on the scenario of formation and evolution of the bulge (e.g. 
Wyse fc Gilmord Il992l ). Mixing of the hot-CNO processed gas forming any GC SG with "pristine" 
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gas is a common re quirement of models i n all clusters show ing the signature of the sodium - oxygen 
anticorrelation (e.g. Prantzos et al. 2007 ; Bekki et al.l 2007 ). If the AGB gas is diluted with a similar 
quantity of accreted metal-rich gas, we require that [Fe/H]~0.44 in the accreting gas. Such a metal 
rich matter itself will have a larg er than solar helium content, say Y ^ 0.30 (Matteucci, private commu- 
nication, see also iRenzinil 119991 ) . so that the very helium rich composition of the mixture forming the 
SG is a tenable hypothesis too. On the contrary, it is difficult to make predictions on abundances of 
the other elements since they will come out from the above described mixture of hot-CNO-processed 
matter (e.g. with low [O/Fe] and large [Na/Fe]) with the bulge iron rich gas. The spe cific abundances 



in the bulge gas depend on its precise evolutio nary history, still not fully understood (IMatteucci et al. 



19991 : iLecureur et all 120071 : iBallero et al||2007j ) 



5. Conclusions 

We have shown that the split HB of Terzan 5 can be interpreted as due to two populations differing 
in helium content and metallicity, and not much different in age (5 (age) lOOMyr). Massive AGB 
stellar models for the chemistry of the FG are compatible with the required helium enhancement, but 
we need that 1) either the AGB matter itself is strongly polluted by SNIa ejecta, before the second 
stage of star formation begins, or 2) the AGB matter is diluted with accreted gas, fastly processed to 
very high metallicity in the bulge stellar environment. This suggestion may help to understand the 
"true" birth of the double population of this cluster, maybe as a mix of age and helium difference in 
the subsequent star formation events. 

We conclude by pointing out that the two alternative scenarios (age or helium difference) predict 
different values for the bright HB and faint HB masses. Specifically while in the merging scenario the 
younger age of the bright HB implies that this population would be ~O.25M more massive than the 
faint HB (in FigJH the mass difference between points C and D), in the scenario proposed in this paper 
the two populations would be almost coeval and their red giant progenitors would have only a small 
mass difference (in the example of Fig. [2J M=O.996M for the bright HB and M=0.979M Q for the faint 
HB). Further dynamical modelling will help to shed further light on the plausibility of the two scenarios 
and on the possible dynamical histories leading to the observed differences in the spatial distribution of 
the two populations. 
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